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Abstract 


Recently, a method was developed of extracting nucleon-nucleon 
(NN) cross sections in the medium directly from experiment. The in- 
medium NN cross sections form the basic ingredients of several heavy- 
ion scattering approaches including the coupled-channel approach 
developed at the Langley Research Center. The ratio of the real to the 
imaginary part of the two-body scattering amplitude in the medium was 
investigated. These ratios are used in combination with the in-medium 
NN cross sections to calculate elastic proton-nucleus cross sections. The 
agreement is excellent with the available experimental data. These cross 
sections are needed for the radiation risk assessment of space missions. 

Introduction 

The transportation of energetic ions in bulk matter is of direct interest in several areas including 
shielding against ions originating from either space radiations or terrestrial accelerators, cosmic ray 
propagation studies in galactic medium, or radiobiological effects resulting from the work place or 
clinical exposures. For carcinogenesis, terrestrial radiation therapy, and radiobiological research, 
knowledge of beam composition and interactions is necessary to properly evaluate the effects on human 
and animal tissues. For the proper assessment of radiation exposures both reliable transport codes and 
accurate input parameters are needed. One such important input is elastic cross sections. The motivation 
of the work is to develop a method for calculating accurate cross sections. These elastic cross sections are 
needed in transport methods both detenninistic and Monte Carlo. 

Nucleon-nucleon (NN) cross sections are the basic ingredients of many approaches (refs. 1 to 10) to 
heavy-ion scattering problem. Most information about these NN cross sections comes from the free two- 
body scattering. These cross sections are significantly modified in a nucleus, due to the presence of other 
nucleons, which is affected through the Pauli exclusion principle and modification of meson field 
coupling constants. (See ref. 11.) Our theoretical approach is based on the coupled-channel method used 
at the Langley Research Center. (See refs. 1 to 6.) This method solves the Schrodinger equation with an 
eikonal approximation. The method needs modifications at low and medium energies. In an earlier work 
(refs. 12 and 13), we developed a unique method of extracting medium modified NN cross sections from 
experiments and found that the renormalization of the free NN cross sections is significant at lower and 
medium energies. These modified in-medium NN cross sections, in combination with the newly 
developed ratio of the real to the imaginary part of the two-body scattering amplitudes in the medium, 
were used to calculate the total cross sections for proton-nucleus collisions (refs. 14 and 15). The blend of 
the renonnalized NN cross sections, the in-medium ratio of the real to the imaginary part of the two-body 
amplitude, and the coupled-channel method gave reliable approach to the total cross sections. The 
purpose of the current paper is threefold: 

1. To put in place a reliable method for calculating elastic cross sections for collisions of protons 
with ions 

2. To use our previously developed NN cross sections in the medium and modified two-body 
ampliUides to calculate elastic cross sections for proton-nucleus collisions 

3. To validate and compare the calculated results with the available experimental data 



4. To provide theoretical results where data are not available (due to nonexistence of experimental 
facilities and/or difficulty in experimental data analysis) 

Method 

For completeness, the essentials of the coupled-channel method are briefly sketched. (See refs. 1 to 6 
for details.) In this approach, the matrix for elastic scattering amplitude is given by 

/(q) = -^ J exp(-Zq.b) {exp[/%(b)]-l} d 2 b (1) 

where 

/and % matrices 

k projectile momentum relative to center of mass 

b projectile impact parameter vector 

q momentum transfer 

X(b) eikonal phase matrix 

The total cross section (J [ot is found from the elastic scattering amplitude by using the optical theorem 
as follows: 

4k 

c tot =y Im[/(<7 = 0 )] (2) 

Equations (1) and (2) give 

a tot = 471 Jo f 1 - ex Pl“ Im (%)] cos [ Re (X)]} b db (3) 

The absorption cross section (CT a [, s ) is given by (refs. 12 and 13) 

poo 

CT abs = 2tt J o (1 - exp[-2Im(x)] } b db (4) 

Having calculated the total and absorption cross sections for many nuclei, elastic cross section is the 
difference of these quantities: 

®el — °tot — ‘Tibs (3) 

Tire eikonal phase matrix x (see refs. 1 to 6 for details) is given by 

X(b) = Xdir(b)-Xex(b) ( 6 ) 
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The direct and exchange tenns are calculated by using the following expressions (refs. 1 to 6): 

Xdk(b) = ^ 21 f exp(tq.b) F^(- q) G«(q) f NN (q) d 2 q (7) 

2k kNN J 

% ex (b) = ApAr f exp(/q.b) F^(-q) G (1) (q) d 2 q 
2k km J 

x — l — f exp(/q'.b) f NN (q + q r ) C(q') d 2 q r (8) 

(2n) z J 

where 

F W and G 1 1 1 projectile and target ground-state one-body fonn factors, respectively 

k NN relative wave number in two-body center-of-mass system 

C correlation function (ref. 6) 

A p and A T mass numbers of projectile and target nuclei, respectively 
The two-body amplitude f NN is parameterized as 

(5(a + i) , B q“ 

/nn ~ — ; k NN ex P — r— (9) 

4tt l 2 

where 

<j two-body cross section 

B slope parameter 

a ratio of real part to imaginary part of forward, two-body amplitude 

It is well-known that the absorption cross section depends on the imaginary part of the eikonal phase 
matrix. This leads us to write the two-body amplitude in the medium f NN m as 

f NN,rn = f rn f NN ( 10 ) 

where f NN is the free NN amplihide and f m is the system- and energy-dependent medium multiplier 
function. (See refs. 12 and 13.) Then the nucleon-nucleon cross sections in the medium ( a NN /fl ) can be 
written as 

C NN,m = fm C NN 
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where a NN is the nucleon-nucleon cross section in free space, and the medium multiplier is given by 


/«=»•! «p(~) + 

where E is the laboratory energy in units of A MeV, D is a parameter in units of MeV, as defined 
subsequently in equation (13). The numbers 12 and 0.14 are in units of MeV and fm -3 , respectively. For 
Aj < 56 (mass number for iron ion representing heavy elements considered in our transport phenomena), 

D = 46.72 + 2.21 A T - (2.25xlO“ 2 )Af (13) 


0.14 


1/3 


exp 


E_ 

D 


( 12 ) 


and for Aj > 57, 


D = 100 MeV 

In equation (12), p av refers to the average density of the colliding system and is 

Pav — A p Pj4t) 

where the density of a nucleus Ay (/ = P,T) is calculated in the hard sphere model and is given by 

p = Ai 

Ai (4tt/3) rf 


(14) 


(15) 


(16) 


where the radius of the nucleus r ; - is defined by 

r\ = 1.29 (r/) mis (17) 

The root-mean- square radius (r ; ) nns is obtained directly from experiment (ref. 16) after “subtraction” of 
the nucleon charge fonn factor (ref. 2). 

Note from equation (3) that total cross section depends on real component of eikonal phase matrix, and 
hence (eqs. (5), (6), and (7)), on the product of aa in two-body amplitude. Since the modification of the 
cross sections in the medium has been determined and tested thoroughly (refs. 12 and 13), the 
modification of a — ratio of real to imaginary part of the two-body amplitude — is studied in the medium 
to calculate the total cross sections (refs. 14 and 15). Some data for total cross sections are available for a 
few systems at high energies. Unfortunately, no data are available for total cross sections at low and 
medium energy range; there are some data for p + Pb in the 100 A MeV range. Therefore, values of the 
medium-modified a have been tested for higher energies. At low and medium energies, our theoretical 
results, which incorporate the in-medium two-body amplitudes, can be validated, if and when 
experimental data become available. 
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A best estimate of medium-modified a takes into account the enhancement of the cross sections 
(ref. 17) and stability and is given by 


a,n 3 exp 


-(£-13A 1/3 ) 2 


K 


5000 


l + exp[(10-E)/75] 


(18) 


where 


K = 0.35 + 0.65 


exp 


— (N-Z) 

3 


(19) 


with N being the neutron number of the nucleus and Z its charge number. 

Equation (3) has also been modified to account for the Coulomb force in the proton-nucleus cross 
sections. This modification has significant effects at low energies and becomes less important as the 
energy increases and practically disappears for energies around 50 A MeV and higher. 


For nucleus-nucleus collisions, the Coulomb energy is given by 


V B 


\MZpZt 

R 


( 20 ) 


where the constant 1.44 is in units of MeV-fm, Zp and Z T are charge numbers for the projectile and target, 
respectively, and R, the radial distance between their centers, is given by 

A 1/3 + A 1/3 

R = r P +r T +l.2- £ — [ j 1 ] — (21) 

Earn 

The number 1.2 in equation (21) is in units of fin- MeV 1/3 . In our earlier work (refs. 12 and 13), these 
expressions were used also for the proton-nucleus collisions to have a unified picture of any colliding 
system. However, as shown in the references, equation (21) overestimates the radial distance between 
proton-nucleus collisions, and hence, equation (20) underestimates the Coulomb energy between them. 
To compensate for this, we multiplied equation (20) by the following factor (refs. 12 and 13), which gives 
the Coulomb multiplier to equation (3): 


( 


Xcoul — 


Cl 




1 + 

V Eon J 


1- 


C 2 V B 


' cm / 


For Ay < 56 (mass number for iron), 


C l = 6.81 - 0. 17A r (1.88 x 10“ 3 )Af 
C 0 = 6.57 - 0.30A r (3.6 x 10" 3 )Af 


( 22 ) 


(23) 
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The constant C\ is in units of MeV. For Aj > 57, 


Ci = 3.0 MeVl 
C 2 = 0.8 J 


(24) 


For the nucleus-nucleus collisions, Cl = 0 MeV and C2 = 1 . This form of Coulomb energy was found to 
work well for the proton-nucleus absorption cross sections (ref. 12). Equation (5) is the main equation. 
The total cross section (eq. (3)) and absorption cross section (eq. (4)) are multiplied by equation (22) to 
get the total and absorption cross sections in the medium and then these are used in equation (5) to get the 
results shown in figures 1 to 6. For clarity, the final expressions used for calculating the total and 
absorption cross sections in the medium are given as 


a to t=4jr| 


1 + Ci 

V E cm J 


1-CoV 


2 V B 


1*00 

J 0 {l-exp[Im(x m )] cos [Re(% w )]} 


bdb 


(25) 


and 


*^abs 271 


/ 1 + C] ^ 

V E cm J 


1-CoV 




\ E cm 


J ’oo 

0 {1 - exp[-2 hn(% m )] ) 


bdb 


(26) 


The constants Cl and C2 are given by equations (23) and (24) and the in-medium eikonal phase matrix 
Xm takes the fonn: 


%m(b) X/w,dir(b) X/w,ex(b) 

The direct and exchange terms in the medium are calculated using the following expressions: 

Ap At 


2Cm.dir(b) 


2n Icnn 


Xm.ex^i 


Ap At 
2k km 


x 


(271) 


J exp(zq.b) C (1) (-q) G (1) (q) /Aw,m(q) d 2 q 

J exp(/q.b) E (1) (-q) G (1) (q) dr q 
— J exp(/q'.b) f N , v>OT (q + q') C(q') d 2 q' 


(27) 


(28) 


(29) 


The two-body amplitude in the medium //v/v,w includes the in-medium nucleon-nucleon cross section (G m ) 


r 0 

JNN,m ~ 

4n 


k NN exp 


' Bi£ ] 
V 2 J 


(30) 


Equations (25) to (30) are used to calculate the elastic cross sections shown in figures 1 to 6. The 
procedure is as follows: Using the in-medium NN cross sections, calculate the two-body amplitude in the 
medium by equation (30). Then use equations (28) and (29) to calculate direct and exchange part of the 
eikonal phase matrix and equation (27) for the full eikonal matrix, which in turn is used in equations (25) 
and (26) to calculate the total and absorption cross sections in the medium. And finally, equation (5) gives 
the results shown in figures 1 to 6. 
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Results and Conclusions 


Figures 1 to 6 show the results of our calculations for the elastic cross sections for proton on 
beryllium, carbon, aluminum, iron, lead, and uranium targets, respectively. The experimental data have 
been taken from the compilation of references 18 and 19. There is paucity of data at lower and 
intermediate energies where the medium modifications play a significant role. For the energy ranges 
considered, where the data are unavailable, our results provide good theoretical values of total cross 
sections because many renormalization effects due to medium, which play an important role in cross 
sections, have been incorporated in the formalism. The reason for the enhancement in elastic cross 
sections in the intermediate energy range is mainly because the real part of the two-body scattering 
amplitude is more dominant compared with the imaginary of the two-body scattering amplitude in this 
energy range. 

We find very good agreement with the experimental results for all the systems at higher energies 
where some data are available. The in-medium cross sections derived earlier in combination with the 
modified ratio of the real to the imaginary part of the amplitude discussed provide good results for the 
proton-nucleus elastic cross sections. It is gratifying to note that the present method gives a consistent 
basic approach for the total reaction and the total cross sections, hence, the total elastic cross sections for 
the entire energy range for all the systems studied. 

The in-medium two-body amplitudes developed in our approach can be used with ease in other nuclear 
processes such as the collision of ions with nuclei as well. The next step is to examine effects on specified 
reaction channels. 
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Figure 2. Elastic cross sections for proton-carbon collision as function of energy. 
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Figure 3. Elastic cross sections for proton-aluminum collision as function of energy. 
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Figure 4. Elastic cross sections for proton-iron collision as function of energy. 
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Figure 6. Elastic cross sections for proton-uranium collision as function of energy. 
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